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ABSTRACT 

The concept of expanding space has come under fire recently as being inadequate 
and even misleading in describing the motion of test particles in the universe. Previous 
investigations have suffered from a number of shortcomings, which we seek to correct. 
We study the motion of test particles in the universe in detail, solving the geodesic 
equations of General Relativity for a number of cosmological models. In particular, 
we use analytic methods to examine whether particles removed from the Hubble flow 
asymptotically rejoin the Hubble flow, a topic that has caused confusion because of 
differing definitions and invalid reasoning. We conclude that particles in eternally 
expanding but otherwise arbitrary universes do not in general rejoin the Hubble flow. 

Key words: cosmology: theory - gravitation - relativity 



1 INTRODUCTION 

The paradigm of modern cosmology is the Friedmann- 
Lema i tre-Robertson- Walker (FLRW) model ( FViedmannl 
Il922l: iLemaitrd Il93ll: iRobertsonI Il935l: IWalkej Il93dl . In 

recent years there has been debate over paradoxical fea- 
tures of the model, and the physical interpretation of 
its dynamics. Attention has been drawn in particular 
to superluminal recession velocities Davis fc Lineweaved 



| 200ll:bavis. Lineweaver. fc Webbll200ll: IChodorowskil2006l: 
ISitnikovl200^ and to the motion of test part i cles in expand- 
ing s pace feavis. Lineweaver. fc We bbll200ll: IWhitindliool 
iGrein fc Elg'ar(2)v....200fi, : .Peacock 2Q06J . It is to the second of 
these issues that we turn our attention. 

The motion of test particles in the FLRW model is 
a fascinating illustration of the interaction between physi- 
cal concepts and quantitative theories. One of the defining 
characteristics of physics is the mathematical precision of 
its predictions. Yet there is more to applying physical laws 
than simply solving equations. In order to make physical 
laws more transparent and accessible, we use physical con- 
cepts that develop an intuition or a mental picture of the 
scenario. A successful physical concept allows us to short- 
cut the mathematics, qualitatively understanding a scenario 
without having to solve the equations. As an example, con- 
sider the north pole of a bar magnet approaching a loop 
of wire. Looking from behind the magnet, we know that 
an anti-clockwise current will be set up in the wire, but 
we need not come to this conclusion by solving Maxwell's 

* E-mail: luke, mfrancis, jbjames, gfl@physics.usyd.edu.au 



equations — Lenz's law will give us the answer without the 
mathematics. 

Attached to the equations of the FLRW model is the 
physical concept that "space is expanding" . Galaxies, we are 
taught, are receding not because they are moving through 
space but because space itself is being stretched between us 
and the galaxy. On the face of it, this concept gives us a good 
intuitive understanding of many cosmological phenomena — 
it helps us understand why the velocity-distance law is linear 
and why light is redshifted as it moves through the uni- 
ve rse. However, it has been attacked r ecently, most notably 
bv lWhitind JgOOi) and lPeacoclj ll200d). as being ina dequate 
in describing local dynamics. Whilst IPeacockl ll2006ll will al- 
low a global form of the expanding space concept — the total 
volume of a closed universe increases with time — he contends 
that: 

there is no local effect on particle dynamics from the global 
expansion of the universe . . . 'Expanding space' is in general a 
dangerously flawed way of thinking about an expanding universe. 

Previous attempts at resolving this debate have suf- 
fered from a number of shortcomings. The most common 
is the overwhelming desire to approximate General Rel- 
ativity (GR) by something else — Newtonian gravity. Spe- 
cial Relativity (SR) or a weak field limit of GR. This is 
probably wrong and certainly unnecessary — why approxi- 
mate when you can use the exact geodesic equations? An- 
other probl em is the small r ange of cosmological models 
considered — iGrOn fc Elgareivl ll2006f) . for example, do a mar- 
vellous job of solving the geodesic equations and then ap- 
ply the solution to just two models, both of which are 
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observationally disfavoured. There is also a dangerous re- 
liance on numerica l calculat i ons t o determine asymptotic 
(t — > oo) behaviour. 'Whiting' ('2004') makes this point about 
iDavis. Lineweaver. fc We bb (200JJ, but doesn't say why an- 
alytic solutions to a Newtonian approximation are better 
than numerical solutions of the exact GR equations. 

In this paper we address these problems by treating 
the asymptotic behaviour of the full GR equations analyti- 
cally, for general (up to non-oscillating, bounded equations 
of state) cosmological models. Section |5| will give the full 
equations of motion for particles in the FLRW model. In 
Section |^ we will analyse test particle motion and correct 
some errors made in previous work. In Section^Jwe consider 
the asymptotic behaviour of particles in the universe, and 
in particular the notion of joining the Hubble flow. Finally, 
in Section |K| we briefly outline the relevance of our results 
for expanding space. Our aim is not to be the flnal word 
in this debate but instead to clarify the relevant features of 
Robertson- Walker spacetime. 



2 THE COSMOLOGICAL GR EQUATIONS OF 
MOTION 

We begin by giving the equations for the FLRW model. The 
Robertson- Walker (RW) metric, which tells us how to mea- 
sure distance and time in a homogeneous and isotropic uni- 
verse, has the line element: 



Sl (de^ -f- sin^ 



where c is the speed of light (hereafter c = 1, unless 
reintroduced for clarity) and ^(^(x) = sinx,XiSinhx for 
k = -fl, 0, — 1, where k is the curvature constant for closed, 
flat and open space respectively. The RW metric relates the 
spacetime interval ds to the cosmic time t and the spheri- 
cal comoving coordinates {x,d,<t>)- More will be said about 
these coordinates later. The scale factor R{t) is the key pre- 
diction of any cosmological model, encapsulating the begin- 
ning, evolution and fate of the universe. We can also deflne 
the Hubble parameter H, which measures the rate of expan- 
sion of the universe: 

In the curved, expanding spacetime of the RW metric, 
we must be very careful when defining distance measures 
(see iLindeil (|l997) for detaUs and Qjpee ( 1Q99) for a sum- 
mary). Throughout this paper, we will use proper distance 
Tp, which is defined as being the radial {d8 = d<j) = 0) space- 
time interval (ds) along a hypersurface of constant cosmic 
time {dt = 0)^. The RW metric then gives the proper dis- 
tance between the origin (x = 0) and x at time t to be: 

^ A thought experiment for measuring proper distance is as fol- 
lows: we imagine being at one end of a giant ruler, pointed at a 
distant object. A volunteer is sent along the ruler to read off the 
distance to the object. Since the universe is expanding, the vol- 
unteer will need to carry a clock that displays cosmic time, and 
note down the time when the measurement was made. When light 
rays have carried the volunteer's result back to us, wc will know 
the proper distance to the object at the time the measurement 
was made. Samuel (2005) criticises proper distance as "violating 
the principle that instantaneous non-local measurements cannot 



rp(t) = R{t)x{t) 



(3) 



The field equations of GR allow us to find R{t) given the 
energy content of the unive rse. The result is the Friedmann 
equations, which, following iHobson. Efstathiou. fc Lasenbvl 
(»2005.i . equation 15.13, we will write as: 



Ro 



-3(1+Wi) 



(4) 



The other two Friedmann eq u ation s can be found in 
iHobson. Efstathiou. fc Lasenbvl (l2005tl . but will not be 
needed here. A subscript zero always refers to a quantity 
evaluated at the present epoch. The sum is over the en- 
ergy components of the universe (labelled i), each with 
corresponding equation of state Wi — pi/pi, where p is 
the pressure and p is the energy density. Equation I^J as- 
sumes that the energy components do not interact and that 
each Wi is a constant^, as it is for most familiar forms 
of energy — matter (w,„ = 0), radiation (wr — 1/3), vac- 
uum energy {wa = —1). The sum includes "curvature en- 
ergy" {ilk,o = -k/(RoHo)^ = 1 - J2i^k^^-o) which has 
Wk = —1/3. This can be thought of as convenient shorthand. 
A particular solution to the Friedmann equations which will 
prove useful in Section 14.21 is the case of a universe with a 
single component w > — 1: 



(5) 



where to is the age of the universe. For w — —1, we have the 
solution R{t) — Ro exp ^^7^^ > where to is the e- folding time 
for the expansion. Wo will not consider phantom energy with 
w < —1 (see Caldwell, Kamionkowski. fc Weinberg 20o3, 
and references therein for details). 

The trajectory of a particle in the universe is computed 
by solving the geodesic equations of GR, given by: 



dV^ dx^dx^^^ 
dA2 dA dA 



(6) 



where A is an affine parameter, the indices a, b, c run from 
to 3 and (x°, ,x^,x^) = {t, x, d, 4>)- The Christoffel symbol 
r^c is given by: 



ra 
hr 



ig'"' {dbQcd + dcQbd 



ddQbc 



For the RW metric, we have 



gab = diag{l,-R ,-R Sk{x),-R Sk{x)sm 9) 



(7) 



(8) 



be made" . This amounts to criticising a spacelike interval for be- 
ing a spacelike interval. In any GR metric, length or distance is 
defined as the spacetime interval along a surface of constant time, 
and as such can never be known instantaneously. This does not 
mean, however, that proper distance is unphysical. It only means 
that it must be reconstructed at a later time from the information 
in light signals. 

■2 Evolving equations of state and interacting components will 
only be considered with regards to their asymptotic behaviour — 
cf. footnote 151 see [Barnes et al. I2005i) and references therein for 
details on such cosmological models. 
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which gives: 

t" + RR (x" + Slix) {9'^ + sin' e<p'^)) =0 
X" + 2Ht'x' - Skix)Ckix) [9'^ + sin' 60") =0 



9 — sin 9 cos t 



+ 2//tV + 2^4^x>' + 2 cot ( 



=0 



^" + sin' 9<i>'^)) 



(9a) 
(9b) 

(9c) 

(9d) 
(9e) 



where a prime means d/dA, and an overdot means A/At. 
Cfe(x) = dS'fc(x)/dx = cosx, l,coshx for k = +1,0,-1. 
Equation H9e^ is the normalisation condition for the four- 
velocity; we set Q = — 1,0, +1 for spacelike, null and time- 
like geodesies respectively. For timelike geodesies, this means 
that the afiine parameter is the proper time measured by a 
clock that follows a given geodesic. We will only consider 
timelike geodesies. 

Equations H9a|l - l|9d(l allow a trivial solution, i.e. Xi 9, 
(f) = constant for all time. Thus, there exists a family of 
free-falling particles that maintain their position in comov- 
ing coordinates — in fact, this is what we mean by comoving 
coordinates. This family of particles defines the Hubble fiow. 

For motion in one dimension, we can choose our 
coordinates to make make the motion purely radial. 
iGrein fc Elgareivl ll200d) study this case and derive the fol- 
lowing useful equations: 



X ={R + CR 

X =Xo ± 



-1/2 

At 



RVl + CR2 



where 



C 



(10a) 
(10b) 

(11) 



is strictly positive. Note that these are equations for x as a 
function of t not A. 



3 MOTION OF TEST PARTICLES 

Following previous work on this topic, we set up the initial 
conditions of our test particle as follows. We place ourselves 
at the origin (x = 0) and the test particle at xo- We con- 
sider the particle to initially have constant proper distance 
rp{to) = i.e. 



At I 







i?(to)xo + i?(to)x(io) = (12) 



The reason for this particular initial condition is quite 
simple. A popular way of visualising expanding space is 
a balloon or a large rubber sheet. Imagine yourself and a 
friend at rest on a large rubber sheet. We cannot directly 
observe spacetime, so we will do this thought experiment in 
the dark^. Suppose you both observe a glowing ball moving 
away from you. "The rubber sheet is being stretched," you 
say. "No it's not," replies your friend, "the sheet is still and 
the ball is rolling away." Together, you come up with an 

^ We will not speculate on how you two came to be standing on 
a rubber sheet in the dark. 



ingenious way of finding out who is right. You take another 
glowing ball, and drop it onto the sheet a certain distance 
away. If the sheet is expanding, then we expect it to carry 
the ball away; if the sheet is still then the recession of the 
first ball was due to a kinematical initial condition. Once 
this is removed, so is the recession. 

The cosmological expansion is a bit more complicated, 
as we have expansion that changes with time due to the self 
gravitation of the energy contents of the universe. We will 
therefore need to consider a range of cosmological models. 
For reasons that will become clear in the next section, we will 
consider models given by Equation 0. We will allow x to 
be negative when a particle passes through the origin, rather 
than have to worry about a change in angular coordinates. 
We have chosen to = 1 and i?o = 1 for each model. We 
started each particle off at comoving coordinate xo = 1/3, 
which ensures that C > 0. Physically, this ensures that we 
do not place the particle beyond the Hubble sphere, which 
would require a velocity relative to the local Hubble fiow 
(peculiar velocity) greater than the speed of light 

The results of solving the geodesic equations are shown 
in Figure The dotted line in the centre panels shows the 
mot ion of the particle as calculated by a Newtonian analysis 
fsee IWhitin'3 ||2004) and lPeacoc3 feOOdl '). The equation of 
motion is: 



Newton n 7") 

r„ = 2Roxo 



1/3 



^oXo 



2/3 



(13) 



The Newtonian result is surprisingly accurate, remaining 
close to the GR solution even up to 100 times the age of the 
universe. It is seen to diverge from the exact solution eventu- 
ally, though, and thus remains only a useful approximation 
for small times. This divergence becomes more apparent as 
we increase xoj i-e. as we approach the Hubble sphere. If 
we consider the w = 0, Einstein de-Sitter universe and set 
Xo = 1, then the solutions d iverge much mo re quickly. Fig- 
ure |5| reproduces Figure 1 of lWhitinS i2004) . overlaying the 
relativistic solution. It is easy to see that whilst the quali- 
tative behaviour is similar, the Newtonian solution is quan- 
titatively different. 

Returning to FigureQ a few points are noteworthy. The 
bottom, leftmost panel {w = —2/3) shows the particle tra- 
jectory moving away from the origin and very quickly be- 
coming indistinguishable from the nearby Hubble flow. The 
other panels do not show the same behaviour, but instead 
the particle moves toward the origin and away on the oppo- 
site side of the sky. Moreover, they don't seem to be attach- 
ing themselves to any particular particle in the Hubble flow. 
But, as noted in the introduction, it is dangerous to try to 
determine asymptotic behaviour from numerical plots. We 
must do it analytically. 



* lGr0n fc Elgaravl l200(J) . equation (22) gives this condition in- 
correctly. The correct formula, which follows directly from their 
equation (21), is xo < 3(1 + «')/2. They then claim to start at 
particle off at xo = 1 in a Milne {w = —1/3) universe, which 
contradicts the previous condition on xo- Therefore, their Figure 
lb) appears to plot a null geodesic. 
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w = -2/3 w = -1/3 w = w=1/3 w = 2/3 




time t 



Figure 1. Comoving radial coordinate and proper distance (solid lines) as a function of time for a particle in radial motion for cosmological 
models with differing values of the equation of state ui given above each column. The dashed lines in the lower panels show the motion 
of nearby particles in the Hubble flow. The dotted line in the centre panels gives the Newtonian solution for the motion of the particle, 
as discussed in the text. Note that the vertical scale changes from panel to panel. 



4 JOINING THE HUBBLE FLOW 



There is disagreement in the hterature as to the fate of 
free particles in an eter naUy expanding universe, i.e . where 
R(t ) ^ 00 as t ^ oo. iDavis. Lineweaver. fc WebbI ll200lll 
and lGron fc Elgareivl (l2006l) claim th at they wiU asy mptot- 
ically rejoin the Hubble flow, whilst IWhitina 1)2004^ states 
that "it cannot be asserted that . . . free particles [are] swept 
into the Hubble flow, even asymptotically." This is an im- 
portant issue because it is often claimed that the expansion 
of space will dampen out all motion through space, so that a 
particle initially removed from the Hubble flow will asymp- 
totically rejoin it. 



4.1 Seven Definitions 

Definition 1 (x ^ Q): A particle with coordinate 
trajectory x(i) asymptotically rejoins the Hubble flow if 
x{t) ^0 as t —> oo. 

The particle is deemed to asymptotically rejoin the Hub- 
ble flow if its velocity through coordinate space approaches 
the velocity through coordinate space of the Hubble flow, 
namely zero. 



On closer inspection, this disagreement stems from dif- 
ferent definitions of what it means to "asymptotically rejoin 
the Hubble flow." In this section we will propose a number 
of precise definitions of this phrase, and then see which ones 
are equivalent and which ones hold in an eternally expanding 
but otherwise arbitrary universe. We will not consider uni- 
verses in which the current expansion becomes a contraction 
at some point in the future. This would be an unnecessary 
and distracting complication when considering the expan- 
sion of space. We do not claim that this list is exhaustive. 



Definition 2 (x ^ Xoo): A particle with coordinate 
trajectory x(i) asymptotically rejoins the Hubble flow if 
X ~* Xoo as t ^ oo, where Xoo is a constant that depends 
on the cosmology and initial conditions. 

The Hubble fiow is defined by having constant coordinates. 
Thus, we consider a particle to approach the Hubble fiow if 
its radial coordinate approaches a constant. The asymptotic 
value of the radial coordinate (Xoo) can be thought of as the 
"rightful place" of the particle in the Hubble fiow. 
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Figure 2. A reproduction of Figure 1 in Whiting 12004), plotting 
proper distance against time for a test particle released from xo = 
1 in an Einstein de-Sitter universe (i.e. flat, matter only). The dot- 
dashed line is the Newtonian solution, whilst the solid line is the 
relativistic solution. The dashed lines represent particles in the 
Hubble flow. The discrepancy between the solutions is obvious. 



Definition 3 (vpec -^Q): A particle with coordinate 
trajectory x{i) asymptotically rejoins the Hubble flow if 
Vpec{t) = R{t)x{t) ^ as t ^ (X. 

We can divide the proper velocity (fp) of a test particle into 
a recession component and a peculiar component as follows: 



(14) 



If we move our coordinate origin so that x(i) = at time 
t, then we see that the proper velocity of the test particle 
is solely its peculiar velocity. Thus peculiar velocity is sim- 
ply proper velocity relative to the local Hubble flow. The 
requirement that Vpc^c(t) — > as t ^ oo is equivalent to the 
velocity relative to the local Hubble flow going to zero. 



Definition 4 frp Vrec): A particle with coordinate 
trajectory asymptotically rejoins the Hubble flow if 
rp — > Vrec{t) as t ^ oo. 

We require that the proper velocity of the test particle ap- 
proaches its recession velocity. This is subtly different from 
definition 3, as will be explained below. To avoid ambigu- 
ity, this definition uses a continuous version of asymptotic 
equivalence: we say that f{x) approaches g{x) as a; ^ oo if 
their ratio approaches unity, i.e. 



/(x) g{x) as a; oo <4> 



1 as x 



(15) 



Definition 5 (Arp ^ 0): A particle with coordinate 
trajectory x(i)) where x{t) — * Xoo as t ^ oo, asymptotically 
rejoins the Hubble flow if Arp = \R{t)Xoc — R(i)x{^)\ ^ 
as f oo. 

Suppose that our test particle is approaching a particular 
coordinate (x = Xoo, cf. definition 2) and that we place a 
reference particle in the Hubble flow at this coordinate. We 



require that the proper distance between the test particle 
and the reference particle approach zero. In other words, 
the test particle sees its rightful place in the Hubble flow 
get closer (in terms of proper distance) asymptotically. Note 
that if we had chosen instead to require that R{t)x{t) 
R(t)Xoa then this definition would have been equivalent to 
Definition 2 by Equation 1151 . 



Definition 6 (zobs —> Zmsm): A particle with observed 
redshift Zobs{tr) at time of reception tr asymptotically 
rejoins the Hubble flow if Zobs{tr) — * Zcosm as tr oo, 
where te is the time of emission of a photon that reaches 
the observer at tr- 

Light emitted from a particle in the Hubble flow is observed 
to be redshifted according to the cosmological redshift for- 
mula: Zcosm = R{tr) / R{te) — 1. For a particle with coordi- 
nate trajectory x(*)j there is an additional Doppler redshift 
resulting from its velocity relative to the Hubble flow: 



1 + Zohs{U) = (1 + 2cosm)(l + ^Dop) 

_ / R{tr)\ ( 1 -I- llpoc(te) 



\R{t,) 



't'pcc e 



(16) 
(17) 



where Upcc is considered positive when the particle's velocity 
through the local Hubble flow points away from us. The 
Doppler term is the familiar redshift of light formula from 
SR, but the formula is derived purely from the RW metric. 

This definition assumes that we are comoving observers 
that measure light signals sent from the test particle. Sup- 
pose at each time we place a reference particle in the Hubble 
flow at the same coordinate as the test particle. The redshift 
of the reference particle, which represents the local Hubble 
flow, will be purely cosmological. This deflnition requires 
that the redshift of the test particle approach the redshift 
of the reference particle. 



Definition 7 (CMB dipole ^ Q): A particle moving 
through a universe containing a cosmic microwave 
background (CMB) asymptotically rejoins the Hubble flow 
if the dipole anisotropy in the CMB goes to zero as t ^ oo. 

In a universe flUed with black-body radiation at a certain 
temperature T, an observer moving through the Hubble flow 
will see the CMB to be hotter in one direction and colder in 
the opposite direction. Indeed, this is exactly what we see 
from Earth — it is known as the "great cosine in the sky" 
and disrupts the isot ropy of the CMB at a level of ~ 10~^ 
iBennett et alj|2003l among others). The maximum differ- 
ence between the temperature as measured by an observer 
in the Hubble flow (To) and ou r test particle with peculiar 
velocity fpec << c is given bv iR^eebles fc Wilkinsonlll968l : 

■ ~3E ' 



iMelchiorri. Melchiorri. fc Signordl2002l) : 



AT 



'-'pec 
C 



(18) 



Finally, note that it is too much to ask that the test 
particle exactly join the Hubble flow after some time, i.e. 
x{t) — Xf, a constant, for all t ^ t{. Only pathological func- 
tions that do not equal their own Taylor series can do this, 
and it is unlikely that such functions will appear as a solu- 
tion to the geodesic equations @. 
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4.2 Comparing the Definitions 

The previous section may appear to be an exercise in pedan- 
tic cosmology, and one iiopes tiiat all the definitions will 
turn out to be equivalent. However, it turns out that only 
three of the above definitions hold in eternally expanding 
but otherwise arbitrary universes, and one of them fails in 
all non-accelerating universes. 

We will need two key results in order to analyse these 
definitions. The first is that all eternally expanding cos- 
mological models approach the single component model of 
Equation as t ^ oo. We can show this directly from 
Equation Q. Consider the universe to contain a number of 
energy components (labelled i), each with constant^ equa- 
tion of state Wi. Now consider the right hand side of Equa- 
tion 0. As t ^ CX3, we know that R{t) — > oo since we 
are only considering eternally expanding universes. Since 
the dependence on R{t) is ggg that, for large 

t, the component with the most negative equation of state 
will dominate the dynamics of the universe. Precisely, let 
the component with the most negative equation of state be 
called the dominant component (i = d), with equation of 
state Wd- Then, for large t 

-3(1+™^) 

(19) 



R{t) 



(20) 



which is Equation Q with w — Wd- For example, in a uni- 
verse with matter density less than critical, we consider 
"curvature energy" to be the dominant component with 
Wd = —1/3. The exact solution for this cosmology (given 
in iHobson. Efstathiou. fc Lasenbvl (|200^, pg. 402) indeed 
shows that R{t) oc t for large t. Thus, Equation H2()|l is a 
general form for R{t) when considering the asymptotic be- 
haviour of the universe^. We can also calculate the deceler- 
ation parameter, q: 



^ = -:RiP = 2(^"'' + ^) 



(21) 



Thus, if Wd > —1/3 then the expansion of universe will 
decelerate; if < —1/3, then the universe will eventually 
accelerate; if Wd ~ —1/3 then the universe will approach a 
coasting universe. 

The second key result is the asymptotic behaviour of 
the integral for x in Equation lllObll when R{t) given by 
Equation 1201 . Whilst the exact indefinite integral unfortu- 
nately involves the hypergeometric function, we can approx- 
imate this function in the limit of large t as by noting that 

^ For a component with an evolving equation of state, consider 
the asymptotic value of the equation of state, i.e. Wi{t) — > Wi^^ as 
t ^ 00. A unbounded equation of state is most likely unphysical. 
An oscillating equation of state will not be considered. 
^ There is a subtlety here — a universe that contains only "cur- 
vature energy" (i.e. an empty universe) is not identical to a uni- 
verse containing a critical density fluid with equation of state 
w = —1/3. Although the dependence of the scale factor on time 
is the same in both universes {R{t) oc t), the empty universe 
has k = —1, whilst the w = —1/3 fluid universe has A; = 0. 
However, Equations jSJ show that if we consider radial geodesies 
(6' = rf)' = 0), then there is no dependence on fc. Thus, the dis- 
tinction between these two universes can be ignored for now. 



Vl + CR^ ~ VCR in this limit. The integral then becomes 
trivial. We now analyse the seven definitions, in order from 
the weakest to the strongest conditions. 

Definition 1 (x ~^ 0): From Equation llObt we can 
see that as t ^ 00 (and R{t) — > 00), x oc R~^ ■ Thus Defi- 
nition 1 holds in all eternally expanding universes, so that 
the velocity of a par ticle through coo r dinat e space will al- 
ways decay to zero. iGrOn fc ElgarOvl (f200ff) use this defi- 
nition when they claim that test particles will rejoin the 
Hubble flow. 



Definition 3 (up 



0): Since Vp^^{t) = R{t)x{t), 



Definition 3 is stronger than Definition 1. However, it still 
holds in all eternally expanding universes, since Equation 
lllObll shows that as f — > 00, VpEc{t) oc R~^ . This is closely 
related to the well known result that the momentum of any 
particle in the universe decays as R~^ . It is this definition 
that is most widely used to justify the claim that test par- 
ticles rejoin the Hubble fiow asymptotically. 

Definition 7 (CMB dipole 0): Equation GHJ shows 
that this definition is equivalent to Definition 3 and thus 
holds in all eternally expanding universes. In particular, this 
shows that the dipole in the CMB will decay faster than the 
CMB temperature itself. 

Definition 2 (x — > Xoo): At first sight. Definition 2 
may appear to be a direct consequence of Definition 1 — 
surely if the derivative of a function approaches zero then 
the function itself will approach a constant. However, it is 
easy to think of counterexamples: f{x) = log(a::). Thus, we 
need to consider Equation lllObll . integrating between to = 1 
and t and using the two approximations discussed at the 
start of this section. This gives: 



dt 



RVT 



CRq 



CR^ 



(22) 
(23) 



where n = ^rn-. r and Ki is a constant that includes xo 

and the primitive of the exact integral evaluated at to- Thus, 
for the integral to be bounded as t ^ 00, we require that: 



2n > 1 



Wd < 1/3 



(24) 



Thus if we use Definition 2 to define what it means to asymp- 
totically rejoin the Hubble flow, then test particles in uni- 
verses where the dominant energy component has equation 
of state Wd ^ 1/3 do not rejoin the Hubble flow. In partic- 
ular, in a universe where the dominant energy component 
is radiation, the comoving coordinate of a test particle re- 
moved from the Hubble fiow increases (or decreases) with- 
out bound. The particle has no rightful place in the Hubble 



flow' 



Definition 4 



Definition 4 appears to be 



identical to Definition 3; surely if rp{t) — Vicc{t)+Vpac{t) and 
Vpoc{t) — > in all eternally expanding universes then Defini- 
tion 4 holds trivially. However, it is possible that Vrac{t) also 



lHartleli200l sets the derivation of Equation llObi as a practice 
problem for undergraduates. The solutions give a very instructive 
derivation using Killing vectors, but then state that the particle 



comes to rest at coordinate x 



given by integrating from 



zero to inflnity. He fails to note that the integral may diverge, so 
that the particle may not come to rest at all. 
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goes to zero, and if it does so as fast or faster than Vpcdt) 
then we are not justified in saying that the proper velocity 
of the test particle approaches its recession velocity. We can 
see this from Equation 1151 : 



rp{t) Vrcc{t) 



Vrcc{t) 

^'rcc(t) + ^'pcc(^) 

Vrcc{t) 
Vpcc{t) Q 



(25) 



(26) 



(27) 



Vroc (t) 

which is stronger than the condition in Definition 3. IWhitin3 
i2004l page 11) hinted at Definition 4: "Peculiar veloci- 
ties do vanish eventually in expanding universes; but so do 
all velocities [emphasis original]", but mistakenly implied 
that it would rule out rejoinin g the Hubble fiow in matter 
domin ated universes {wd = Ol. lDavis. Lineweaver. fc WebbI 
i200ll) use Definition 4 (see their equation (11) and follow- 
ing) but mistakenly believe that it follows automatically 
from the success of Definition 3. 

As we noted previously, ast^oo, Vpccit) oc R^^ . Thus: 

Vp^c{t) (X R~'^ (X t~" (28) 

We now consider «roc(i) = R{t)x{t). Using Equation 1231 we 
approximate x by: 



X{t) ^K^± 



if n / 1/2 



VCR^ [logt if n= 1/2 
T^,"^^^T^ ifrj/1/2 

Vrcc(t) OC K2t + K:i < 

U-^/^logt ifn = l/2 



(29) 
(30) 



where the Ki will be used keep track of constants. Now, we 
need to consider three cases: 

• If n > 1/2, then n — 1 > —n . Thus the dominant term 
in Equation 13011 is the first term so that Viec{t) oc for 
large t. Then, as t oo: 



Vrec{t) 







(31) 



Thus when n > 1/2, Upcc(i) goes to zero faster than i'rcc(t), 
meaning that Definition 4 holds. 

• If n = 1/2, then Vrcc{t) oc t~^''^ logt. Then, as t ^ oo: 

,-1/2 

(32) 



Vpec{t) 



Vrec{t) t-l/2 1ogt 

Thus Definition 4 holds when n = 1/2. 

• If n < 1/2, then n — 1 < — n so that Vrcc(t) oc t~". 
Then, as t ^ oo: 



Vpcc{t) t 
_L L_L (2<i — 

Vrcc{t) 



(33) 



Thus when n < 1/2, Vrec{t) and i'rcc(i) approach zero at the 
same rate. It is not true that rp{t) — > Vicc{t) as t — > oo and 
Definition 4 fails in this case. 

If we use Definition 4 to define what it means to asymp- 
totically rejoin the Hubble fiow, then test particles in uni- 
verses where the dominant energy component has equation 
of state Wd > 1/3 do not rejoin the Hubble fiow. Note that 
Definition 4 holds in a universe where the dominant energy 
component is radiation, unlike Definition 2, which fails. 



Definition 6 (zobs — > Zcosm): Once again we are 
tempted to assume that the success of Definition 3 will en- 
sure that this definition will hold in all universes. The argu- 
ment proceeds as before: if 



1 -I- ZDop = 



1 + ^'pcc(te 
1 ' '^pcc (^e 



(34) 



and Vpcc{t) ^ in all eternally expanding universes then 
Definition 6 holds trivially. And our caveat is the same — we 
must be careful of the case where Zcosm also goes to zero. 

To do this, we need to express Zobs in terms of reception 
time tr and then consider the limit tr ^ oo. In fact it is much 
easier to express everything in terms of the emission time 
and then consider te — > oo; te < tr guarantees that both 
cases will have identical limiting behaviour. 

In the limit of small Upoc, we have that 

2^Dop ^ ^pcc OC 

using Equation 1281 . Also, 1 -I- Zcosm = = -jk, so 

that the task at hand is to express tr in terms of te- Consider 
a light ray travelling along a null (ds — 0), ingoing {dx < 0), 
radial {d6 = d(j> — 0) geodesic. From the RW metric: 



Xite) - Xitr) = 



dt = - R{t)dx 
dt 



R{t) 



(35) 
(36) 



where we now place the receiver at the origin: xi^r) = 0. For 
the case of a test particle moved removed from the Hub- 
ble fiow, x(^e) is given by Equation llObL approximated 
by Equation 1291 . An immediate consequence of combining 
these equations is that Definition 6 must work in acceler- 
ating and coasting universes {wd ^ —1/3), since in these 
universes Zcosm does not go to zero. Thus we need only con- 
sider < n < 1. We will leave the n = 1/2 case to the 
reader: it involves log t as with Definition 4. 

With this in mind, we can derive an expression for tr'. 



tr 



which leads to the following expression for 1 + Zc. 

1 + Zcosm = (l + Kit-'' + A-sie"''""'"* ' 

^1 + K^tr + ifrtr''""' 



(37) 

(38) 
(39) 



where the last expression is calculated using the binomial 
theorem for non- integer exponents*. This leads to the fol- 
lowing three cases: 

• If n > 1/2, then n— 1 > — n so that Zcosm oc Then, 
as t — > oo: 



2Don t 

oc 



Zcosm t^ 







(40) 



Thus when n > 1/2, zoop goes to zero faster than Zcosm, 
meaning that Definition 6 holds. 

• If = 1/2, it turns out that Zcosm oc t~^^^ logt. Then, 
as t oo: 



ZDon t 

— oc 



-1/2 







Zcosm t-^/'^logt 

Thus Definition 4 holds when n = 1/2. 
* See IWeissteirJ l2006l) for a reminder. 



(41) 
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Figure 3. Redshift as a function of time for cosmological models with w = (left) and w = I (right) (examples of fluids with ui = 1 are free 
massless scalar fields and shear energy, such as superhorizon gravitational waves; see lLind'eJ ll997l) 'l. The dashed line is the cosmological 
redshift (zcosm), the dotted line is the Doppler redshift (zoop) a^nd the solid line is the observed redshift (zobs = — l+(l+2cosm)(l+2Dop))- 
The left panel shows the asymptotic dominance of the cosmological redshift for to 1/3, whilst the right panel shows the identical 
asymptotic behaviour of the cosmological and Doppler redshift at large t for w > 1/3. 



• If n < 1/2, then n — 1 < — n so that Zc. 
as t ^ oo: 



oc t 



~ 

^cosm t 



1 



Then, 



(42) 



Thus when n < 1/2, Zcosm and zoop approach zero at the 
same rate. It is not true that Zobs — > .Zcosm and Definition 6 
fails in this case. 

If we use Definition 6 to define what it means to asymp- 
totically rejoin the Hubble flow, then test particles in uni- 
verses where the dominant energy component has equation 
of state Wd > 1/3 do not rejoin the Hubble flow. The situa- 
tion is illustrated in Figure|H] When w = the Doppler red- 
shift decays away much faster than the cosmological redshift, 
so that Definition 6 holds. However, when lo = 1 the Doppler 
and the cosmological redshift decay at the same rate; which 
one is greater depends on the cosmological model and the 
initial conditions (hidden in the constants Ki). Thus, Defi- 
nition 6 fails in this universe. 

Definition 5 (Arp — > 0); We know already that this 
definition fails in some cases since it relies on Definition 2. 
Thus we begin with the assumption that Wd < 1/3, i.e. 
n > 1/2. Now, we have that: 



Arp =\R{t)xo. - R{t)x{t)\ 
dt 



--R{t) 

-Kit) 



f 



At 



At 



to Ry/TTCW 



oc t 



rVTTcW 

At_ 

f2n 



(43) 
(44) 

(45) 

(46) 
(47) 



Hence the requirement that Arp ^ as t —> oo is only 
met for n > 1, i.e. Wd < —1/3 or g < 0. In particular, for 
a universe where the dominant component has Wd = —1/3 , 
Arp approaches a constant. For example. IWhitin3 i2004ll . 



page 10 reaches this conclusion for an underdense, matter 
only universe. If we use Definition 5 to define what it means 
to asymptotically rejoin the Hubble flow, then test parti- 
cles in universes where the dominant energy component has 
equation of state ^ — 1/3 do not rejoin the Hubble flow. 
This means that particles asymptotically rejoin the Hubble 
flow only in universes that eventually accelerate. Definition 
5 is the stronge st of the definitions, and is the one used by 
IWhitin j l)2004h . page 10, to justify the claim that particles 
in a matter dominated universe {wd = 0) do not join the 
Hubble fiow. 

The failure of this definition can be illustrated by look- 
ing again at Figure^ The first column {w — —2/3) clearly 
shows a particle that satisfies Definition 5. However, in the 
second and third columns {w — —1/3, 0), in the bottom pan- 
els, the lowest dashed line is the trajectory of the particle 
in the Hubble flow whose comoving coordinate is equal to 
Xoo, i.e. the rightful place of the particle in the Hubble flow. 
In the centre column, the test particle trajectory and the 
rightful place trajectory move apart, whilst in the second 
column they are separated by a constant. Thus, the parti- 
cle never joins the Hubble flow in the sense of its trajectory 
being indistinguishable from the trajectory of a particle in 
the Hubble flow. 

A summary of the different deflnitions and the condi- 
tions for their fulfilment is given in Tabled 



5 PARTICLE MOTION, THE HUBBLE FLOW 
AND EXPANDING SPACE 

The attack on the physical concept of expanding space has 
centred on the motion of test particles in the universe, as 
discussed in Section |S1 Whiting (2004), Peacock (2006) and 
others claim expanding space fails to adequately explain test 
particle motion because we expect that expanding space will 
carry the particle away from the origin, as stretching rubber 
would carry the glowing ball of Section |3 away. Whilst we 
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Table 1. Seven definitions and the conditions for their fulfilment. 





Definition 


Hold in all 
expanding 
universes? 


If no, conditions 

on Wd for 
definition to hold 


1 




yes 




2 




no 


Wd < 1/3 


3 


Vpcc 


yes 




4 


^ ^rcc 


no 


Wd ^ 1/3 


5 


Arp -* 


no 


Wd < -1/3 


6 


^ohs * -^cosm 


no 


Wd ^ 1/3 


7 


CMB dipole 


yes 





will not attempt to resolve this debate here, we believe that 
this may be a misunderstanding of expanding space that is 
fostered by a flaw in the rubber sheet analogy. Particles in 
the Hubble flow do not feel any force, as they are free-falling. 
Thus, thinking of expanding space as a frictional or viscous 
force (like objects on a rubber sheet) is incorrect. 

Finally, how are we to understand the failure of so many 
of the deflnitions of joining the Hubble flow in Section 
We contend that expanding space predicts t hat the peculia r 
velocity of a test particle will approach zero jPeacoc kl200eD : 

through the accumulated Lorentz transforms required to over- 
take successively more distant particles that are moving with the 
Hubble flow. 

Thus we expect Deflnition 3 (along with Definitions 3 and 7, 
which are weaker or equivalent conditions) to hold in all eter- 
nally expanding universes. However, expanding space does 
not lead us to expect that Definitions 2, 4, 5 and 6 will hold. 
We contend that the problem is not that expanding space 
has mislead us, but that describing the decay of Vpcc as join- 
ing the Hubble flow is a misnomer. The correspondence be- 
tween this particular characteristic of the trajectory of the 
test particle and the trajectory of particles in the Hubble 
flow leads us to expect that all of the features of the trajec- 
tory (X) X) ^p> '^pi Vpcc) will approach those of the Hubble 
flow trajectories (respectively, x ~ Xoo, x = 0, r-p = R{t)xoo, 
i'p = Vrcc, Vpcc = 0). But this is too much to ask from the 
expansion of the universe. As we have seen, many of these 
conditions depend on the acceleration and deceleration of 
the universe, rather than just its expansion. 

We contend that the correct definition of asymptoti- 
cally rejoining the Hubble flow is that all the features of the 
test particle trajectory approach the corresponding features 
of the Hubble flow trajectories. Selecting one feature of the 
trajectory on which to base our definition is arbitrary and 
leads to a multitude of conflicting claims. All seven defini- 
tions have equal claim to the title of the definition of joining 
the Hubble flow. We therefore require that all deflnitions 
hold, which is equivalent to just requiring Deflnition 5 to 
hold, as it is the strongest definition. It follows that it is not 
a general feature of expanding universes that test particles 
asymptotically rejoin the Hubble flow. 
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